1/1 


AD-A139  138  EXTREME  VALUE  THEORV  AND  DEPENDENCE (U>  NORTH  CAROLINA 
UNIV  AT  CHAPEL  HILL  CENTER  FOR  STOCHASTIC  PROCESSES 
H  R  LEHDBETTER  JUL  85  TR-189  AF0SR-TR-85-8692 
UNCLASSIFIED  F49628-82-C-8888  F/G  12/1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  - 1963 ~ A 


AD-A159 


CENTER  FOR  STOCHASTIC  PROCESSES 


Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina 


Extreme  value  theory  and  dependence 


M.R.  Leadbetter 


.  'sr 


Technical  Report  No.  109 
July  1985 


DT1C 

fcEUECTE 
X  SEP  1  3  1985 


9L 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

lb.'  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  OlSTRIBUTION/AVAILABILlTY  of  REPORT 

UNLIMITEE^Proved  for  public  release  ; 
distribution  unlimited. 

2b.  OECLASSIFICATION/OOWNGRAOING  SCHEOULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER  IS) 

Technical  Rept.  NO.  109 

5.  MONITORING  ORGANIZATION  REPORT  NUMBERiS) 

SFOSR-TR-  S  3-  0  6  92 

6a.  NAME  OF  “SRPORMING  ORGANIZATION  |6o.  QFPiCS  SYMBOL 

iff  applicable r 

Center  for  Stochastic  Processes 

7a.  NAME  OF  MONITORING  ORGANIZATION 

! 

6c.  A.OORESS  <Cify.  State  ana  ZIP  Code/ 

Statistics  Dept.,  321  PH  039A 

UNC 

Chapel  Hill,  NC  27514 

|  7b.  aOORESS  City.  State  and  ZIP  Coae> 

]  afosr/nm 

2££fM»«0  «>332-64« 

a».  name  of  funoing/sponsoring 

ORGANIZATION 

AFOSR 

Bb.  OFFICE  SYMBOL 
if  f  applicable ' 

8c.  AOORESS  <Cit y.  State  and  ZIP  Codet 

Bolling  AFB 

Washington,  DC  20332 

11.  *!7La  tfnciuae  Security  Classification) 

EXTREME  VALUE  THEORY  AND 

DEPENDENCE 

|  to.  SOURCE  OF  funding  NCS. 

PROGRAM 

I  PROJECT  ! 

TAS< 

WORK  UNIT 

element  no. 

NO. 

NO. 

NO. 

1 

J _ 

i  i 

A5 

1 

l 

1 

12.  PERSONAL  AUTHOR(S) 


13a  TYPE  of  REPORT 

technical 


is.  Supplementary  notation 


M.R.  Leadbetter 


13b.  tim£  COVEREO 

eaoM  9/84 _ 


1«.  DATE  OF  REPORT  iYr..  Mo..  Oayj  15.  PAGE  COUNT 

July  1985  6 


17.  COSATI  COOES 


FIELD  GROUP  1  SUB.  GR. 


18.  SUBJECT  TERMS  'Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Stochastic  processes,  extreme  value  theory,  high  level 
exceedances 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  &y  bloc *  number) 

—  The  purpose  of  this  paper  is  to  give  a  very  brief  account  of  some  of  the  essential  ideas 
underlying  classical  extreme  value  theory,  and  to  see  how  these  are  used  (modified  as  neces¬ 
sary)  for  dependent  cases.  In  particular  it  will  be  shown  how  the  classical  theory  still 
applies  for^moderately  dependentl,*‘s'tat ionary  sequences,  but  that  under  higher  local  depen¬ 
dence,  clustering  of  high  values  occurs,  requiring  modifications  of  the  theory  especially 
as  it  involves  order  statistics  other  than  the  maximum.  Underlying  concepts  (.especially 
point  process  convergence  results)  are  emphasized.  «... fQ_  L _  ■ 


'Jn  'Y  ' 1 


Z jiy  i  v  a  \  v*.  f /  j. 


20.  CiSTRi  BUT, ON.  A  v  Ai  lA3i  l.  TY  OF  ABSTRACT 

-N;_„C3.=  ED'LNl  mi'ED  1?  SAME  as  RP’-  _  3T,C  .SE  RS  _ 


,21  ABSTRACT  SECURITY  CLASSIFICATION 

'  UNCLASSIFIED 


name  of  respcns  blE  nCivixjuai. 


J/Ua 


EXTREME  VALUE  THEORY  AND  DEPENDENCE 


by  M.R.  Leadbetter 


K^'rv'frion  For 

NT  IS  GRA&I 
DTIC  TAB 
Ur announced 
J  ■  Jtif 1 cat  ion _ 


T4 


□ 


By-  - _  _ 

Distribution/ 

_ Availability  Codes 

(Avail  and/or 
3ist  |  Special 


Oepartiaent  of  Statistics  i  | 

University  of  North  Carolina  /j -'A j 

Chapel  Hill,  N.C.  27514  fl  j 

Abstract.  The  purpose  of  this  paper  is  to  give  a  very  brief  account  of  some  of  tne  1 
essential  ideas  underlying  classical  extreme  value  theory,  and  to  see  how  these  are 
used  (modified  as  necessary)  for  dependent  cases.  In  particular  it  will  be  shown  how 
the  classical  theory  still  applies  for  "moderately  dependent"  stationary  sequences, 
but  that  under  higher  local  dependence,  clustering  of  high  values  occurs,  requiring 
modifications  of  tt#e  theory  especially  as  it  involves  order  statistics  other  than  the 
maximum.  Underlying  concepts  (especially  point  process  convergence  results)  are 
emphasized. 


COPY 
.  INSPtCTf 


1.  Notation,  and  the  Classical  Theory. 

Throughout  we  write  Mn  ■  max(E1,C2« • • • *or 

any  sequence  (£n>  of  random  variables.  If  the 

^  are  i.i.d.  with  marginal  d.f.  F  then  Mr  has 

the  d.f.  P(M  <x)«Fn(x)  so  that  the  distribut- 
n  — 

ion  of  Mn  is  precisely  known  if  F  is  known. 
However  in  practice  F  is  not  necessarily  known 
precisely,  and  approximations  less  dependent  on 
the  exact  form  of  F  are  useful.  In  this  vein  a 
central  contribution  of  the  classical  theory  is 
the  following  result  which  restricts  the  "type" 
of  limiting  d.f.  which  may  have  (under  linear 
normalizations) : 


Theorem  1,1  (Extremal  Types  Theorem).  Let 
Mn = max(Ci,E2, . . . ,Cn) .  where  are  i.i.d. 
random  variables.  If  for  some  constants 


a  >  O.b  ,  we  have 
n  ’  n’ 

(1.1)  P(an(Mn  -  b  )  <  x}  *  G(x) 
n  n  n  — 

for  some  non-degenerate  d.f.  G,  then  G  is  one 
of  the  "three  extreme  value  types": 

Type  I:  G(x)  *  exp(-e~*),  -o»<x<«>; 

Type  II:  G(x)  *  .  -cu  ,  _ 

(family)  \exp(-x  ) ,  for  some  a >  0, 


according  as  x  <  0  or  x  >  0 


Type  III:  G(x)  .  j«P(- (-*)“>.  f°r  some  cr>0, 

(family)  according  as  x  <_  0  or  x  >  0 
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In  these  d.f.'s  x,may  be  replaced  by  ax  ♦  b  for 
any  a>0,b.  This  is  consistent  with  the  term- 
iology  we  use  that  two  d.f.'s  G,H  are  of  the 
same  type  if  G(x)  ■  H(ax+b)  for  some  a>0,b. 

The  classical  theory  provides  domain  of  attract¬ 
ion  criteria  determining  the  type  of  limit  in 
terms  of  the  general  behaviour  of  the  tail  1-F 
of  the  d.f .  F  (see  [5])  .  Thus  Equation  (1.1)  may 
often  be  used  to  provide  an  approximation  for 

the  distribution  of  M  when  n  is  large.  Obvious- 
n 

ly  relevant  questions  of  rates  of  convergence 
and  estimation  of  the  normalizing  constants  in 
practice  are  not  part  of  this  paper  (see  [S] 
and  references  therein  for  discussions  of 
these  topics) . 

While  no  detailed  proof  of  Theorem  1.1  can  be 
given  here,  it  is  useful  to  see  the  principal 
ideas  involved.  These  can  be  grouped  into  the 
following  two  packages: 

Package  I:  The  class  of  d.f.'s  G  which  may 
occur  as  limits  in  Theorem  1.1  are  precisely 

the  max  stable  d.f.'s  i.e.  d.f.’s  G  such  that 

k  k 

G  is  of  the  same  type  as  G  (G  (x)  * G(a^x  ♦  6^) 

for  some  >  0,6) ,  for  each  k«l,2 . 

Package  II:  A  d.f.  G  is  max  stable  if  and  only 
if  it  is  one  of  the  three  extreme  value  types. 

The  result  of  Package  I  follows  trivially  for 
the  i.i.d.  sequence  considered  by  noting  that 


.*  . 


• 


(1.1)  may  be  written  as  Fn(x/a  +b  )-*-G(x)  from 

n  n 

which  it  follows,  on  replacing  n  by  n k  ,  that 

(1.2)  P{ank(Mn  -  bnfc)  £  x)  ■*  G1^k(x)  as  n  +  « 

for  each  k*  1,2,...  .  A  well  known  result  of 

Khintchine  (cf.  (5  Theorem  1 . 2. 3]) states  that  if 

a  sequence  (here  )  has  limits  under  more 

than  one  1  inear  normalization,  they  must  be  of 

1/k 

the  same  type,  so  that  G  is  of  the  same  type 
as  G,  from  which  max  stability  is  immediate. 

Package  II  requires  a  longer  proof  but  it 
should  be  noted  that  this  is  totally  non-probab- 
ilistic  and  does  not  involve  the  original  i.i.d. 
sequence  in  any  way  (see  (5  Theorem  1.4.1]). 
Hence  Package  II  applies  equally  well  to  non 
i.i.d.  situations. 

Theorem  1.1  is  the  central  result  of  the 
classical  theory,  though  the  use  of  linear 
normalizations  is  coming  more  into  question. 

In  fact  (1.1)  may  be  written  as 

(1.3)  P{M  <  u  }-*-e'T  as  n-*-“°. 

n  —  n 

where  u  =  x/a  *b  ,  x  «  -logG(x).  This  suggests 
n  n  n 

the  possibility  of  using  e.g.  functions  un * 
un(x)  which  are  not  necessarily  linear  in  x. 
Further,  for  any  sequence  of  constants  {un> 
(dependent  on  a  parameter  x  or  not)  (1.3)  may 
be  written  as  Fn(un)-*e"T  and  by  taking  logs 
one  sees  simply  that 

(1.4)  n [  1  -  F (un) ]  •*  t  as  n-*«®  . 

The  converse  result  is  just  as  easily  shown 
giving  the  following  simple  but  important 
lemma . 

Lemma  1 . 2  Let  (un)  be  any  real  sequence,  and 
0  i  T  1  ”•  Then  (1.3)  and  (1.4)  are  equivalent 
for  the  i.i.d.  sequence  (£n)  with  d.f.  F. 

We  shall  use  the  notation  u^fx)  to  denote  any 

sequence  satisfying  (1.4).  It  is  of  interest 

(and  makes  the  notation  more  compelling)  to 

note  that  if  u  (t)  exists  for  one  t  with 
n 

0  <  x  <  °°  ,  it  exists  for  all  such  x  (and  this 

happens  if  and  only  if  (1  -  F(x))/(1  -  F(x-))  -*■  1 

as  x-*-x  ($ee  [S,  Theorem  1.7. 13]) ,  where 
F 


xp(£  ®)  is  the  upper  endpoint  of  the  d.f.  F.) 

Lemma  1.2  is  central  in  the  consideration  of 

domains  of  attraction.  But  it  is  also  basic  in 

showing  how  the  asymptotic  distributions  of 

extreme  order  statistics  are  determined  by  that 

for  the  maximum.  To  see  this  let  S  denote  the 

n 

number  of  exceedances  of  the  level  u  by 

Cj,...,Cn*  the  number  of  i  such  that 

^  >un>  1  £  i  ±  n.  Sn  is  clearly  a  binomial 

r.v.  with  parameters  (n,  p„ *  1  -  F (uR) )  so  that 

if  (1.4)  holds,  np  x  and  S  has  a  Poisson 
n  n 

limit  with  mean  x,  i.e.  P{S  ■  r)  -*■  e*Txr/r  ! 

n 

fkl 

Now  it  is  easily  seen  that  if  M  is  the  kth 

largest  of  ^^  ....  £n>  then  M(h) when  there 

are  no  more  than  k-1  values  of  £.  greater  than 

u_,  i.e.  S  <k-l,  so  that 
n  n  — - 

k-1 

P{M*k)  <  un)  *  P{S  <  k-1)  -  l  e'V/r! 

r»0 

In  particular  if  (1.1)  holds  (so  that  has 
the  asymptotic  distribution  G)  then  the  ident¬ 
ification  u  *  x/a  +b  x  =  -log  G (x)  and  Lemma 
n  n  n 

1.2  show  that 
‘  m 

(1.5)  Pta^M^  -  bn)  <  x)-G(x)  l  (-log G(x)) r/rl 

r«0 

Thus  the  asymptotic  distribution  for  M  deter- 
(k)  n 

mines  that  for  each  Mv  (using  the  same 
n 

normalizing  constants  as  Mn) . 

The  Poisson  property  of  exceedances  may  be 

pursued  further  with  advantage.  Specifically 

let  u^  ■  un(x)  satisfy  (1.4)  for  some  x >  0. 

Write  Nn  for  the  point  process  on  the  interval 

[0,1]  having  an  atom  (event)  at  each  point  i/n 

for  which  That  is  a  time  scale  change 

by  1/n  is  made,  and  converts  an  exceedance  at 

i  (Ct  >  u„)  to  a  point  at  i/n  in  [0,1].  may 

be  called  the  exceedance  point  process  for  the 

level  u  •  If  N  (B)  denotes  the  number  of  atoms 
n  n 

of  Nn  in  the  subset  B  of  [0,1],  and  (1.4)  holds, 
it  is  readily  shown  that  if  B  is  an  interval 
then  N  (B)  is  asymptotically  Poisson  with  mean 
X  r,(B)  where  m(B)  is  the  length  of  B.  Further 
N^(Bj)  and  are  clearly  independent  when 

B^  and  B2  are  disjoint.  This  suggests  that  the 


point  process  is  talcing  on  a  Poisson  character 
as  n  increases  and  indeed  it  may  be  shown  by 
such  consideration  that  if  (1.4)  holds  then  full 
weak  convergence  holds, 
d 

(1.6)  Nn  ♦  N  as  n-*  ® 

where  N  is  a  Poisson  Process  on  [0,1]  with 
intensity  t. 

A  harvest  of  corollaries  may  be  reaped  from  this 
and  related  results  involving  more  than  one 
level.  In  particular  the  asymptotic  joint  dis¬ 
tribution  of  any  group  of  extreme  order  statis¬ 
tics  and  their  locations  may  be  obtained.  In 
the  i.i.d.  case  these  results  may  be  also 
obtained  directly  but  the  point  process  approach 
is  illuminating,  and  most  useful  in  dependent 
cases  to  be  considered  next. 


2.  Stationary  Sequences  and  the  Extremal  Types 
Theorem. 

Turning  now  to  dependent  situations  we  consider 
the  case  of  a  which  is  strictly  stationary 

in  the  standard  sense  that  its  finite  dimension¬ 
al  distributions  F .  (x....x  )■ 

ll  'p  1  P 

P{tAi  <  Xi---Cip  <  xp}  «e  such  that  F^.  ...yj  = 

F.  :  for  any  choice  of  p,  i,...i  ,  j. 

V  p  1  P 

Other  cases  (e.g.  Markov-cf.  [7],  [8])  may  be 
treated,  but  stationary  sequences  are  adequate 
to  illustrate  the  effects  of  dependence. 

As  far  as  the  Extremal  Types  Theorem  is  concer¬ 
ned,  Loynes  ([6])  took  the  first  and  most 
significant  step  away  from  independence  by 
showing  that  the  result  remains  true  under 
strong  mixing  assumptions.  That  is  if  {£n}  is 
stationary  and  strongly  mixing  and  (1.1)  holds 
for  some  non-degenerate  G,  then  G  must  be  of 
extreme  value  type.  It  is  obvious  from  Loynes' 
proof  that  the  full  force  of  strong  mixing  is 
not  required  and  weaker  assumptions  will  suffice. 
The  following  condition  "D(un)"  is  defined  with 
respect  to  a  given  sequence  {u  )  and  is  con¬ 
venient  and  useful.  (It  is  clearly  possible  to 
weaken  D(un)  very  slightly  for  the  present 


purpose  -  as  has  been  made  explicit  in  [7],  but 

D(u  )  is  useful  in  other  contexts  also.  In  the 
n 

following  statement  F . (u)  is  written  for 

11  *P 

...j.  (u...u). 

1  P 

Definition.  If  {u^}  is  any  real  sequence,  the 

condition  D(un)  is  said  to  hold  if  for  any 

l<i,...<i  <j,<j,...<  j  ,<n,j,-i  >  l  we 
-  1  -  p  J2  p  —  1  p - 

have  |F.  (u  ) 

1  p  1 1  V 

-  F  .  (u  )F.  .  (u  ) |  <  a  where 

if  •  •  •  i  n  1 . « •  •  i  i  n  n ,  i 

1  p  Jp* 

a  ,  -*•  0  for  some  sequence  £  *  o(n) . 
it,  t  n 

n 

D(u^)  mildly  restricts  long  range  dependence 
(and  is  satisfied  e.g.  for  stationary  normal 
sequences  under  the  very  weak  covariance  con¬ 
dition  rn  log  n-»  0.  A  detailed  discussion  of 
the  condition  may  be  found  in  [5]). 

The  Extremal  Types  Theorem  then  says  that  if 

{£n>  is  strictly  stationary  and  satisfies  (1.1) 

and  if  D(un)  holds  for  all  sequences  of  the 

form  u  *  x/a  *b  (-co<x<°°)  then  G  is  again  of 
n  n  n 

extreme  value  type.  The  method  of  proof  is 
simply  to  show  that  such  a  G  must  be  max  stable 
and  hence  of  extreme  value  type  as  in  the  i.i.d. 
case.  The  proof  of  max  stability  rests  on  the 
following  basic  result: 

Lemma  2.1  If  D(u^)  holds  for  the  stationary 
sequence  {£n)  (and  a  given  real  sequence  {u^}) 
then  P{M  <  u  }  -  P^{Mr  .  c  u  }  -►  0  for  any 
k-  1,2,...  . 

It  follows  simply  from  (1.1)  by  using  this 

lemma  (with  nk  replacing  n  and  identifying  un 

with  x/a  +b  )  that  (1.2)  holds  and  hence  G  is 
n  n 

max  stable  exactly  as  in  the  i.i.d.  case.  The 
proof  of  the  lemma  is  achieved  by  a  standard 
type  of  argument  used  to  "reduce  dependent  to 
independent  cases  (cf.  [6]  and  earlier  papers 
in  dependent  central  limit  theory).  Very 
roughly  in  this  case  the  integers  l...n  are 
divided  into  k  consecutive  groups  and  approximate 
independence  of  the  maxima  on  each  is  used. 

This  approximate  independence  is  established  via 
D(u^)  by  "snipping"  an  expanding  but  relatively 


of  e. 


saall  piece  froa  each  group  to  give  the  separa¬ 
tion  required  for  D(un),  (cf.  [  5  Section  3.2] 
for  details). 


3.  The  Effect  of  Dependence  on  the  Asymptotic 

Distribution  of  N  . 

n 

Again  let  be  stationary  with  marginal  d.f. 

F.  Following  Loynes,  define  the  associated 

independent  sequence  to  be  an  i.i.d.  sequence 

{£  )  with  the  same  d.f.  F.  Write  M  * 
n  ^  n 

max(£. • • •  £_) •  It  is  natural  to  ask  whether 
has  an  asymptotic  distribution  if  Mn  does, 
and  convesely  and,  if  both  do,  whether  they  are 
of  the  same  extremal  type  or  related  in  some 
specific  way.  This  matter  can  be  resolved  in 
unexpectedly  explicit  ways  for  most  stationary 
sequences,  in  terms  of  a  single  parameter. 

As  before  let  u  (x)  satisfy  (1.4)  for  each 

5  -r 

t  >  0  so  that  P{M  <  u  (t))  +  e  by  Lemma  1.2. 
n  —  n 

Now  Loynes  ([6])  showed  that  for  strongly 
mixing  sequences,  if  lim  P{M  £u  (x) )  exists 
for  all  x  it  must  have  the  form  e'  for  some 
9,  0  <  0  <  1.  This  may  also  be  shown  under 
DCu^)  -  assumptions  and  indeed  if  the  limit 

exists  for  one  x  it  exists  for  all  x,  and  is 

*  0T 

e  for  some  9.  This  0  may  be  called  the 
extremal  index  of  the  process  and  exists  under 
wide  conditions.  For  i.i.d.  (and  many  station¬ 
ary)  sequences  9*1,  but  all  values  of  9  in 
[0,1]  are  possible,  though  9*0  is  rather 
pathological.  9*1  for  a  stationary  normal 
process  satisfying  the  covariance  condition 
r^  log  n-*0.  An  example  of  a  case  with  0<9<  1 
is  the  following. 

Example  3. 1  Let  be  i.i.d.  with  d.f.  H 

and  write  f, .  *  max(n  .  ,n  .  ,).  Then  (5  )  is 
j  J  j*l  n 

stationary  with  d.f.  F*H2  and  easy  calculation 

shows  that  if  u  (x)  satisfies  (1.4)  then 
n 

n [ 1  -  H(u^(r)) |  *  x/2  and 
PIMn-Un(T)'  =  P  max(ll...nn)  iun(T)}P{Vl 
iun(x)}-e*T/2 

so  that  (F  )  has  extremal  index  9*  1/2. 
n 

Modifications  of  this  example  yield  other  values 


The  following  result  from  [4]  illustrates  the 
use  of  the  index  0. 

Theorem  3. 2  Suppose  that  the  stationary 

sequence  (£n)  has  extremal  index  0 >  0.  Let 

(vn)  be  any  constants  and  0<_P£l.  Then 

P{fi  <  v  }  -*■  p  if  and  only  if  P(M  <  v  }  p0  . 
n  —  n  n  —  n 

This  clearly  exhibits  the  relations  between  M 
^  n 
and  M  .  For  example  if  there  is  an  "i.i.d. 

limit"  G  P{a  (M  -  b  )  <  x}-»G(x),  then  we  can 
n  n  n  —  ' 

write  v  =  x/a  +  b  and  obtain  P{a  (M  -  b  )  <  x} 
Qn  nn  n  nn- 

-*■  G°(x).  Thus  the  limiting  d.f.  for  M  is 

0  ** 
just  G  .  This  is  of  the  same  type  as  G 

Q 

(G  (x)*G(ax  +  b))  as  is  easily  verified  for 
each  extremal  distribution.  The  converse  also 
holds  so  that  if  9>0,Mn  has  an  asymptotic 
distribution  if  and  only  if  does,  they  have 
the  same  normalizing  constants, and  the  limits 
are  of  the  same  type  ( indeed  powers  of  each 
other).  If  9  *  1  the  limits  are  identical  (see 
[4]  for  a  complete  discussion). 


4.  Poisson 'Results  and  Extreme  Order  Statistics 
when  0*1. 

When  9*1,  the  Poisson  results  for  exceedances 
go  through  as  in  the  i.i.d.  case.  In  parti¬ 
cular,  the  number  S  of  exceedances  of  u 
n  n 

satisfying  (1.4)  by  C ^ - -Cn  is  asymptotically 

Poisson,  leading  to  the  same  asymptotic  distri- 
(k) 

but  ion  (1.5)  for  M  *  as  in  the  i.i.d.  case, 
n 

This  Poisson  limit  is  best  shown  as  a  corollary 
of  Poisson  convergence  of  the  exceedance  point 
process  Nn  defined  on  [0,1]  as  in  Section  1, 
with  points  at  those  i/n  for  which  >  %• 

Specifically  if  un*un(x)  satisfies  (1.4), 

P{M  <  u  }  -*•  e"T  and  hence  P{M  <  u  }  e~T  by 
n  —  n  n  —  n 

Theorem  3.2  (0=1).  But  the  events  [M  <u  ), 

n  —  n 

Nn ( ( 0 , 1  ] )  =0  are  clearly  equivalent  so  that 
P{Nn([0, 1] )  »  0}  ■*  e  T.  The  same  argument  in 
fact  shows  easily  that  P(Nn(B)  *  0)  •*  e‘Tn*W 
where  m  is  Lebesgue  measure,  for  any  finite 
union  B  of  subintervals  of  [0,1],  so  that 
P{Nn(B)  *  0) -*•  P{N(B)  =  0}  for  such  sets  B,  where 
N  is  a  Poisson  Process  on  [0,1]  with  intensity 


•  •  v  v  v  s*  s'  -. 


-  ■-.*.•  ■w- 


V-'.  • 


t.  This  and  the  easily  proved  fact  that 

E  N  (B)  -*  tu(B)  «  E  N(B)  are  sufficient  to  show 
n  d 

that  N  -  N  in  the  full  weak  convergence 
n 

sense,  by  a  theorem  of  Kallenberg  (3,  Theorem 
4.7). 

In  particular  ((0,1]))  is  asymptotically 

Poisson  with  mean  t.  If  has  the  asymptotic 

distribution  G  as  in  (1.1)  it  follows  readily 

by  the  usual  identification  u  * x/a  ♦  b 

n  n  n  ... 

t  =  -  log  G(x)  that  the  kth  largest  values  Mn 

again  have  the  classical  forms  (l.S).  Indeed, 

it  may  be  shown  by  considering  exceedances  of 

more  than  one  level  that  joint  distributions  of 
(k) 

the  (and  their  locations)  have  their  i.i.d. 

n 

forms  when  0=1. 

(k) 

5.  Clustering  and  its  effect  on  when 

0<  0  <  1. 

The  above  discussion  shows  that  when  0=1  the 
classical  results  hold  without  change.  The 
less  usual  -  but  still  non-pathological  -  case 
is  where  0<  0<  1,  which  we  now  consider.  In 
this  case,  which  typically  involves  higher 
local  dependence  of  the  sequence  terms,  clust¬ 
ering  of  exceedances  may  occur.  For  instance, 
it  is  readily  seen  that  in  Example  3.1,  exceed¬ 
ances  of  un(r)  by  occur  in  (at  least) 
pairs.  Correspondingly  it  is  found  that  the 
limiting  point  process  has  double  points, 
occurring  at  positions  which  form  a  Poisson 
Process  with  intensity  t/2. 


(2) 

However  e.g.  Mn  may  be  the  second  largest 

value  in  the  cluster  when  M  occurs,  or  the 

n 

largest  in  some  other  cluster,  so  that  cluster 
structure  becomes  important. 

In  fact  if  the  cluster  size  in  the  limiting 
Compound  Poisson  Process  has  distribution 
ff(i),  i*  1,2...  and  if  has  the  asymptotic 
distribution  G  as  in  (1.1),  then  the  asymptotic 
distribution  for  the  kth  largest  is 

'irf-v;1' 

-GOOll.Y  Y  <-l°»  <H«»J  „•)(!)) 
j*l  i=j  jl 

*  1 

where  it  J(i)  is  the  j-fold  convolution  of  the 
cluster  size  distribution  ir(i).  (These  latter 
may  be  writter  as  limits  of  probabilities 
associated  with  the  C-^s). 

For  completeness  we  state  the  main  theorem 
(due  to  Hsing)  on  which  this  is  based.  This 
requires  a  modest  strengthening  of  the  D(un) 
condition  to  a  form  A(un)  which  will  not  be 
specified  here  (but  is  still  of  similar  type, 
and  widely  applicable  -  see  [1]  or  [2]  for  the 
details  and  proof  of  the  theorem). 

Theorem  5.1  Suppose  t  >  0  is  constant  and 
A(un(x))  holds.  If  the  exceedance  point  process 
Nn  converges  in  distribution  to  some  point 

process  N,  the  later  must  be  Compound  Poisson 
with  Laplace  transform 

-/fdN  -ex/J  fl-L(f(t))]dt 

E  e  *  e 


In  general  when  0<  0<  1  exceedances  tend  to 
occur  in  groups  or  clusters  of  stochastic  size 
and  the  limiting  point  process  is  "Compound 
Poisson"  consisting  of  multiple  events  of 
random  size  occurring  at  points  of  [0, 1J  which 
form  a  Poisson  Process  with  intensity  t/2. 

As  noted  already  the  only  effect  on  Mn  of 
taking  (0  <)0<  1  is  to  replace  the  "i.i.d. 
limit  G"  by  G®.  However,  the  distribution  of 
other  order  statistics  are  more  radically 

altered.  This  may  be  understood  intuitively 
by  noting  that  the  maximum  is  simply  the 
maximum  of  the  greatest  values  in  each  cluster. 


where  L  is  the  Laplace  Transform  of  some  pro¬ 
bability  distribution  n  on  (1,2,...)  and 
0< 0 <  1  . 

Throughout  this  paper  we  have  considered  point 
processes  of  exceedance  of  one  level.  However 
as  noted,  joint  distributions  of  order  statis¬ 
tics  may  be  discussed  by  consideration  of 
"vector  valued"  point  processes  involving 
exceedances  of  e.g.  k  levels  u  (t,)  ...  u  (t)  . 
Alternatively  a  so  called  "complete  convergence 
result"  may  be  used  under  slightly  more 
restrictive  conditions)  to  obtain  these  results. 
Such  results  typically  involve  point  processes 


in  the  plane  consisting  of  points  at  (j/n, 

u*l(£.))  where  u'1  is  the  inverse  function  of 
n  j  n 

u^(t).  In  simple  cases  a  Poisson  limit  holds 
(in  the  plane)  leading  to  all  the  joint 
asymptotic  distributions  of  order  statistics 
(cf.  (S,  Sec.  5.7],  In  cases  with  aore  local 
dependence  (0  <  1)  the  liait  involves  a 
point  process  with  clustering  which  can  be  ex¬ 
plicitly  defined  cf.((l])in  a  manner  analogous 
to  Theorea  5.1. 
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